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Abstract 

The origin of the highest-energy cosmic rays remains elusive. The decay of a 
superheavy particle (X) into an ultra-energetic neutrino which scatters from 
a relic (anti-)neutrino at the Z-resonance has attractive features. Given the 
necessary X mass of 10 14 ^ 15 GeV, the required lifetime, 10 15 ~ 16 y, renders 
model-building a serious challenge but three logical possibilities are consid- 
ered: (i) X is a Higgs scalar in SU (15) belonging to high-rank representa- 
tion, leading to power-enhanced lifetime; (ii) a global X quantum number has 
exponentially-suppressed symmetry-breaking by instantons; and (iii) with ad- 
ditional space dimension(s) localisation of X within the real-world brane leads 
to gaussian decay suppression, the most efficient of the suppression mecha- 
nisms considered. 
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The confluence of cosmology and particle phenomenology benefits both disciplines and 
can lead to important new insights. 

For protons propagating through the cosmological background radiation there is an en- 
ergy cut-off, as discussed in e.g. ]1], well-known as the GKZ effect [0|§, at an energy of 
E ~ 5 x 10 19 eV. Above this energy, the photoproduction of pions at the 3-3 resonance 
provides an energy attenuation that prohibits travel over a distance greater than ~ 50 Mpc. 

Nevertheless, air showers initiated by a proton (or photon) with energies above the GKZ 
bound have been observed @], and this fact needs explanation. 

One possibility is that the origin involves the decay of a superheavy particle as in |J (for 
related earlier works, see e.g. [|]||) but that the decay now produces a high energy neutrino 
which scatters from a relic background neutrino at the Z pole and produces the primary. 
This Z-burst scenario was suggested in [0] and further analysed in . 

The kinematics of the neutrino-neutrino collision at the Z pole requires an energy 
Eresonance = M (Zf / (2M (v)) and taking M(u) = 0.07 eV as suggested by the Su- 
perKamiokande data gives M resonance ~ 10 23 eV, just as needed to explain the data. This is 
the most attractive feature of the model. 

We first estimate the mass and lifetime of the superheavy particle needed to fit the data. 
This requires two relationships derived in ||. Namely the flux of cosmic rays beyond the 
GKZ cut-off is estimated as: 
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For the numerical dimensionless coefficients we find the values C\ = 0.33 and C2 = 12.8 
which we use in the following analysis. The notation is: iV is the number of protons and 



photons per annihilation event; i] is the relic neutrino density relative to the present photon 
number density r\ = (n V)re u c /nry )0 ); Qx is the contribution of X particles to the energy 
density, relative to the critical density; h is the Hubble constant in units of 100km/ s/Mpc; 
B v is the branching ratio of X into neutrinos; to is the age of the universe; tx is the lifetime 
of X; M(v) is the neutrino mass; and M(X) is the mass of X. 

Assuming central values of all other parameters we plot the allowed region of M(X) and 
Tx in Figure 1; variations in N,rj,Qx, h, B u , Miv) can extend the allowed region but here 
we need only the order of magnitude estimate. 

The value of M(X) must certainly exceed 2E resonance so that a two body decay can lead 
to a neutrino acquiring enough energy. Higher energies can be red-shifted down to E resonance 
if the progenitor X particle is at a red shift z > 0. 

The resultant spectrum will cut-off at M(X)/2 and will be expected to provide a two- 
component type of overall spectrum, with a dip around E ~ Eqkz a s can be seen in the 
data @. 

Since Z decay gives rise to approximately 10 times as many photons as nucleons the model 
predicts a concomitant number of high-energy photons as cosmic-ray primaries. Because the 



data is sparse, it is not yet possible to discriminate on this basis, as discussed in |10|; this is 
an important prediction of the Z-burst scenario. 

From the above analysis we conclude that the required particle properties M(X) and 
Tx for the hypothetical state X are well defined in order of magnitude. Namely, the mass 
M(X) should lie between 10 14 and 10 15 GeV and the lifetime t(X) should lie between 10 16 
and 10 17 years. 

The remainder of the paper will discuss three possible microscopic theories or, better, 
scenarios for this combination of M(X) and t(X). We present these three scenarios in what 
we regard as their increasing appeal, from (i) to (iii). 



(i) Power suppression. 

The expectation for a particle of this mass is that, unless it is absolutely stable due 
to some exact conservation law, it will decay exceedingly quickly with a lifetime expected 
to be r < 10 -24 seconds. Since the required lifetime is larger by some 46 or so orders of 
magnitude, the longevity is the principal difficulty, as emphasized in ]IJ. 

One viewpoint is that this extraordinary suppression of the decay is an argument against 
the model, as is the problem, already mentioned, of super-high-energy photons concomitant 
with the protons. 



Let us here take the viewpoint, as discussed in |T(J that the photons are a prediction of 
the model, to be tested in future experiments, rather than a fatal flaw. The data on HECR 
is probably too sparse to reach any stronger conclusion. 

Therefore the only remaining question is longevity. 

The first scenario (i) is that considered (in a different model) in 0. We assume the 
particle X is a boson and posit a coupling 

^fZ:::Z ^ ) ) ) 

where the power is p = 3(n — l)/2. Let us assume that such a coupling is gravity-induced 
and that M is the reduced Planck mass M ~ M(Pl) ~ 10 18 GeV. Then one expects the 
lifetime t(X) to be of the order of magnitude 

in which the mass ratio is jjfjpfj — — 10~ 4 . To arrive at a suppression of 10~ 46 thus 
requires 2p ~ 12 — 15 and n ~ 5 — 6. Thus the X field must have a high tensorial rank. If 
this is too much for the reader, skip to scenario (ii). 

In PJ the case n = 2 was considered. In the spontaneous breaking of S77(15) theory 



IIJ such a tensor appears "naturally" in the Higgs sector. There is no apparent need for 
such a high rank as n = 5 or 6, but equally no reason for their absence. The dimensions of 
such scalar representations in SU(15) are astronomical - even for n=2 the dimension JT|.p~2 



is 14,175 while for n=5 and 6 this becomes respectively 125,846,784 and 1,367,127,216. 



It is difficult to believe that such power suppression could be responsible for the longevity. 
It is a logical possibility which appears highly contrived. Thus exponential or gaussian 
suppression is more appealing. 



(ii) Exponential suppression. 

Let us assume that the superheavy particle X carries a conserved quantum number Qx 
(analagous to baryon number, B) and that in perturbation theory the quantum number is 
exactly conserved. If there is no open channel which conserves Qx then the state will be 
absolutely stable. 



In the case of B in the standard model, it was first shown in 1976 by 't Hooft ||T3|JT4 
that nonperturbative instanton effects violate conservation and lead to decay of otherwise 
stable states such as the proton. The resutant rate is typically exponentially suppressed 
by an exponential of the form exp(— constant/ g 2 ) where g is the gauge coupling constant. 



Many other examples of such suppression are covered in [jnj . 

Thus, one scenario is that Qx generates a symmetry of the lagrangian but X decays with 
exponential suppression due to instanton effects. We mention this only for completeness - 
any quantitative estimation would require many hypotheses. 

(iii) Gaussian suppression. 

This scenario which is, in our opinion, the most appealing involves the assumption of at 
least one extra spatial dimension. We will take five space-time dimensions, four space and 
one time. 

Let the coordinates be (xo, x, y) with y as the hypothetical extra dimension on which we 
now focus. 

It used to be thought, up to a decade ago, that any such y must be compactified at 
or beyond the GUT scale of (10 l6 GeV)- 1 ~ KT 32 m (recall l(GeV)- 1 ~ 2x 10 16 m). 
In 1990, Antoniadis JTB| was the first to entertain very much larger compactification scales 



^1 TeV) 1 ~ 10 19 m. In 1998 it was pointed out ||T7|-pO|] that, although the strong 



and electroweak interactions of the quarks and leptons need be confined to a region of y not 
exceeding 10~ 19 m (the real- world brane on which we live), the gravitational interaction could 
be decompactified even out to 1 mm = 10 -3 m without contradicting experimental data, 
offering the possibility of detecting such additional dimensions by deviation from Newton's 
Law of Gravity at millimeter scales. 

Models in which the fifth dimension contains a real-world brane and a suitably sepa- 
rated Planck brane which delimits gravitational propagation have been discussed in pl| , P2| . 
Such models are of interest mainly because they suggest how to incorporate gravity in the 
conformality approach [f23| -|28H which ab initio describes a flat (gravitationless) space-time. 

Let us assume, therefore, that the standard model states are all confined within a real- 
world brane with a thickness of order 10~ 19 m in the y direction. Following |29| (see also 
SOfl and, for the many fold universe, |31[) a scalar field ^(x^,y) which has a $ domain wall 



in the y dimension. In the vicinity of the wall centered by convention at y = 0, and with 



the normalizations of |29|] the field has the value 

*(y) = 2fi 2 y (3) 

First consider chiral fermions (after all, X could be a fermion but we will consider the 
boson possibility later). In this case we write the five-dimensional action 

S = J d 4 xdy^i[i^d/dy^ + $>{y) - mj*i + ... (4) 

The fermion is now localised at y = m,/ (2fi 2 ). If the Higgs H(y) is unlocalised inside the 
domain wall then the resulting coupling of to and H has a gaussian suppression 

exp(—C (mi — mj) 2 ) / fx 2 ) (5) 

where C is a coefficient of order unity. This is the gaussian overlap of the gaussian tails of the 
two wave functions. The thickness of the real- world brane is r ~ (a0 _1 while the separation 
of the two fermion wave functions is a ~ (Ay)ij. To obtain the required suppression of 
~ 10~ 46 we need r/cr ~ 10. For example, if r ~ 10 -19 m, one needs o ~ 10 -20 m. Clearly 
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only the ratio r/a matters, but need not be a large number, in order to obtain the necessary 
suppression. 

It is worth remarking that this gaussian suppression of Yukawa couplings by localization 
in the fifth dimension has a long historical counterpart in the localization of states on orbifold 



The superheavy particle X may not be a fermion but a boson, e.g. the Higgs scalar 
considered in scenario (i) above. Fortunately we can easily extend the localization argument 
of |29| to the case of a boson and obtain a similar result for the gaussian suppression of the 
decay amplitude and consequent longevity. We replace Eq.(§) by the following: 



and, analagously to the fermion ^, the boson <fr is localized around y = ms/Zfi 2 when $ 
again has a domain wall centered at y — 0. Identifying cf) with X then provides the required 
longevity by the same mechanism. It would be amusing if the highest-energy cosmic rays 
provided the first evidence for an extra spatial dimension! 

To summarize the Z-burst mechanism for the highest energy cosmic rays, it has two 
positive features: 

1) The resonant energy derived from the Z mass and the SuperKamiokande neutrino 
mass is numerically close to the required energy. 

2) The spectrum is predicted to have the two-component shape suggested by the present 
data. 

On the other hand, there are also two apparently negative features: 

3) The concomitant high-energy photons are not confirmed by present data. Better data 
will confirm or refute this important prediction. 

4) The longevity of the superheavy particle is such a challenge to microscpic model- 
building that it may render the model less credible. 

It is point (4) which we have attempted to ameliorate in the present article. 



singularities starting with [|3"2}-|3"4|] in 1987 and subsequent derivative literature [jHSj-fRJ . 




(6) 
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Figure 1. 

Allowed region of M(X) — Tx from Eq.(^) and Eq.(H) of the text. In the Figure ax = 
Tx(10 7 t )~ 1 where t is the age of the universe and bx = M(X)/(10 u GeV). Variations in 
N, r], fix, h, B u , M{y) can extend the allowed region but we use only such order of magnitude 
estimates. 
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